Fractional dynamics from the ordinary Langevin equation 
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We consider the usual Langevin equation depending on an internal time. This parameter is 
substituted by a first passage time of a self-similar Markov process. Then the Gaussian process is 
parent, and the hitting time process is directing. The probability to find the resulting process at 
the real time is defined by the integral relationship between the probability densities of the parent 
and directing processes. The corresponding master equation becomes the fractional Fokker-Planck 
equation. We show that the resulting process has non-Markovian properties, all its moments are 
finite, the fluctuation-dissipation relation and the H-theorem hold. 

PACS numbers: 05.40.-a, 05.60.-k, 05.40.Fb 



I. INTRODUCTION 

The Langevin equation is a powerful tool for the study 
of dynamical properties of many interesting systems in 
physics, chemistry and engineering [l|, The success 
of the approach rests on the description of macroscopic 
quantities starting from microscopic dynamics, where the 
effect of fast degrees of freedom (heat bath in statistical 
physics and solid state physics, short wavelength modes 
in meteorological and climate models, etc) can be often 
taken into account by noise 0-Q . Gaussian noise leads to 
normal diffusion with a mean square displacement that 
grows linearly in time and to an exponential relaxation. 
Equivalently, the phenomena can be also described by 
the ordinary Fokker-Planck equation (FPE) for the time 
evolution of the probability density of the random pro- 
cesses. 

However, many systems exhibit anomalous behavior in 
their transport and relaxation properties 0, Q • Anoma- 
lous diffusion has the mean square displacement increas- 
ing as a (non-linear) power law in time, and anomalous 
relaxation shows a slow power law decay in the long-time 
limit. The attempt to state a dynamical foundation in 
statistical physics, as well as the great interest in un- 
derstanding the physical mechanism leading to anoma- 
lous diffusion/relaxation, calls into being the generalized 
Langevin equations [1, ■ The generalizations affect ei- 
ther the equation form itself (for example, via the mem- 
ory kernel) or/and the character of correlations in the 
fluctuating force The way to the description of 

anomalous diffusion/relaxation is not unique. In this pa- 
per we show that the ordinary Langevin equation can re- 
sult in anomalous diffusion/relaxation owing to the fact 
that the temporal degree of freedom becomes stochastic. 
The approach clarifies a microscopic derivation and inter- 
pretation for the fractional FPE. The ordinary Langevin 
equation is a particular case of the new model. 

The paper is organized as follows. In Sec. |ll]we intro- 
duce the concept of the stochastic time clock. The new 



clock (random process) generalizes the deterministic time 
clock of the ordinary Langevin equation and governs by 
the random process described by the stochastic differen- 
tial equation. The directing process arises from a self- 
similar a-stable random process of temporal steps. Us- 
ing properties of the stochastic time clock, in Sec. IIIII we 
write the corresponding master equation with the frac- 
tional derivative of time. To know the solution of the 
usual FPE with a time-independent kernel, one can find 
immediately the solution of its fractional generalization 
via the integral relation. The ordinary Langevin equation 
has a stationary state. The same feature remain valid, 
when passing to the stochastic time clock. For this case. 
Sec. llVl is devoted to the fluctuation-dissipation relation 
and the H-theorem. The randomization of time clock can 
be also applicable for the general kinetic equation (Sec. 
W\ . The fact is illustrated on a concrete example, the 
relaxation in a two-state system. We end the paper with 
a short summary in Sec. I VII 



II. STOCHASTIC TIME ARROW 

The main feature of time is its direction. Time is only 
running from the past to the future. In our consider- 
ation we intend to save the property of time. For the 
ordinary Langevin equation the time variable is deter- 
ministic. Now set this variable as an internal parameter 
T. The motion of a point particle of velocity V(t) in a 
thermal bath is determined by a viscous friction 7V and 
random collisions W{t), by means of 



dV{T) = ~-iV{t) dT + dWir) 



(1) 
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As usual, W(t) is a Wiener process with zero mean and 
variance per unit of time equal to 2D. Let us random- 
ize the time clock of the process V{t). No every ran- 
dom process is suitable for our goal. First of all the 
appropriate process must be strictly non-decreasing. As- 
sume that the time variable is a sum of random tempo- 
ral intervals T^. Let Ti be independent identically dis- 
tributed variables. It is not necessary to know the exact 
form of their probability distribution. Their belonging 



2 



10 



Time, sec 



10 



FIG. 1: Single realization of a continuous time random walk 
with random waiting-times Ti between successive random 
jumps of velocity Vi. 



to the strict domain of attraction of a a-stable distri- 
bution (0 < q; < 1) is quite enough. The parameter 
restriction < a < 1 arises from the need to keep the 
random time steps Ti as non-negative random variables. 
The sum of random variables n~^/"(Ti + • • ■ + T„), n e N 
converges in distribution to the a-stable one. As has 
been shown in , there exists the limit of the following 
process T^^(t) = {[r/ArJ Elli"^"^^' under 

Ar 0, where r is the internal time separated on dis- 
crete values with a step Ar, and [aj denotes the integer 
part of a. Let us now make use of the limit passage from 
"discrete steps" to "continuous ones". The new process 
satisfies the relation T(t) f t'^/°'T{1), where ^ means 
equal in distribution. The position vector of a walking 
particle at the true time t is defined by the number of 
jumps up to time t. This discrete counting process is 
Nt = max{n G N | — 0- The continuous limit 

of the discrete counting process {Nt]t>Q is the hitting 
time process ^(t) = inf{a; | T{x) > The hit- 

ting time S{t) is called also a first passage time. For a 
fixed time it represents the first passage of the stochas- 
tic time evolution above this time level. The random 
process S{t) is just non-decreasing and depends on the 
true time t. We choose it for a new time clock (stochastic 
time arrow), assuming its statistical independence on the 
random variable V. 

Although the random process S{t) is self-similar, it 
has neither stationary nor independent increments, and 
all its moments are finite [3 [lH . This process is non- 
Markovian, but it is inverse to the continuous limit of 
a Markov random process of temporal steps r(r), i. e. 
S{T{t)) — T. The analytical form of the probability 
density of the random variable S{t) can be calculated as 



follows. According to [1^, the expectation (e 



-vS{t) 



Eai—vt"). After the Laplace transform of the Mittag- 
LefHer function with respect to t, the expectation can 
be easily inverted analytically with respect v. Then the 
probability density of the process S{t) is written as 



(2) 



Br 



where Br denotes the Bromwich path. This probabil- 
ity density characterizes the probability to be at the 
internal time t on the real time t. After the variable 
transform ut ^ u and denoting z = r/t"', the func- 
tion p^{t,T) takes the form t~"Fa{z), where Fa{z) = 
2^ /g^ gU-zu deforming the Bromwich 

path into the Hankel path, we find the Taylor series of 
the function Fa{z), i. e. 
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k=0 



k\T{l 



ka) 



(3) 



where r(a;) is the usual Gamma function. Since the ra- 
dius of convergence of the power series (jS]) can be proven 
to be infinite for < a < 1, the function Fa{z) is entire 
in z. Thus, the exchange between the series and the in- 
tegral in the calculations of the Taylor series ^ is quite 
legitimate. The Laplace image of the function Fa{z) is 
expressed in terms of the Mittag-Leffler function 



Fa{z) dz = Ea{-C) , Z > 0. 



x) is equal to the Mittag-Leffler function 



Feller conjectured and Pollard proved in 1948 that the 
Mittag-Leffler function £;Q(-a;) = E,?=o(~^)"/r(l+"") 
is completely monotonic for x > 0, if < a < 1. More- 
over, Ea(—x) is an entire function of order 1/a for a > 
pjif . Hence, by Feller [l^, one can conclude that the 
function Fa(z) is non- negative in z > 0. Taking into ac- 
count the normalization relation Fa{z) dz = 1, the 
function p^{t, r) is just a probability density. If a = 1/2, 
from the series ([3]) it is easy to recognize the well-known 
Gaussian function Fa{z) = tt^^/^ exp(— z^/4). It should 
be observed here that the basic Cauchy and Signaling 
problems of the time fractional diffusion-wave equation 
can be expressed in terms of the function Fa{z) (see more 
details in [13). 



III. FRACTIONAL FOKKER-PLANCK 
EQUATION 

The Langevin equation directed by the stochastic time 
clock S{t) can be written in the form 

dV{S{t)) = ~-fV{S{t)) dS{t) + dW{S{t)). (4) 

The statement is quite justified. The process V{S{t)) is 
a continuous martingale, and the directing process S{t) is 
a continuous submartingale with respect to an appropri- 
ate filtration [l^. All the moments of both parent and 
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directing processes are finite. The same concerns the 
process V{S{t)). Now the random walk of a particle is 
defined by two Markov processes, random waiting times 
T(t) between random jumps V^(t). The discrete exam- 
ple of such a walk is shown in Fig.[TJ A similar approach, 
i. e. the modeling of anomalous diffusion by two indepen- 
dent random processes indexed by a common continuous 
parameter, has been already suggested in 18]. However, 
the inverse process to the time evolution was not com- 
pletely defined. It was not recognized as a first passage 
process. No direct relationship between Eqs. ([T|) and (|3]) 
was found. This is also a main difference between the 
subject of ll9|] and our paper. 

The resulting process V{S{t)) is subordinated to V{t), 
called the parent process, and is directed by S{t), called 
the directing process The hitting time process S{t) 
just satisfies necessary properties imposed on any direct- 
ing process (independent, positive and non-decreasing). 
The directing process is often referred to as the random- 
ized time or operational time. In other words, the sub- 
ordinated process V{S{t)) is obtained by randomizing 
the time clock of a random process VJt) using a random 
process S{t). In compliance with [l^j, the probability 
density of the random process Vt — V{S{t)) is expressed 
by the integral relation 



(5) 



where p^ (V, r) represents the probability to find the ran- 
dom value V on the internal time r. Recall that p^{t, r) 
is the probability to be at the operational time r on the 
real time t. It is well known that the stochastic differen- 
tial equation of type ([T]) is equivalent to the correspond- 
ing FPE. In particular, the probability density p^{V,t) 
obeys the standard FPE 



dp^{V,T) 
dr 



Lfpp^(V,t), 



where Lpp is a time-independent Fokker-Planck oper- 
ator whose exact form is not important here. The 
Laplace transform of the function p^*{V,t) with respect 
to time replaces the integral relation ([5]) by the algebraic 



one, p"'(y,w) 



,a-l-V 



p {VjU"). Acting the operator 



Lpp on the Laplace image, we obtain [ipp p"'](V, u) = 
P"" (y, u) -QiV) 'u"-\ where Q{V) is the initial con- 
dition. The inverse Laplace transform gives the fractional 
FPE 

p^^{V,t)^ Q{V) + 



1 



dr(t-r)"-MiFPP"'](V,T). (6) 



r(a) JO 

This equation has also the equivalent form 
9>^'(y,t) Q{V)t 



r(i-a) 



LFPp^^{V,t)., 



where d°'/dt°' denotes the Liouville-Riemann fractional 
differential operator of order a [lO). Our analysis gen- 
eralizes the mathematical treatment of [l^ and shows 
that the Fokker-Planck operator can have a more general 
form rather than only with a temporally constant force 
field. Another approach to the description of anoma- 
lous transport in external fields is developed in [2l'| . The 
consideration is based on a generalization of the classical 
Chapman-Kolmogorov equation. An interesting justifica- 
tion of the generalized Chapman-Kolmogorov equation is 
that trapping events are superimposed on the Langevin 
dynamics, with a waiting time distribution with infinite 
mean. By the choice of special forms for the transfer ker- 
nel and the probability density function of the waiting 
time between any two successive jump events in the gen- 
eralized equation, one can recover some models discussed 
in the literature. 

If the probability density p^ (V, r) is known explicitly, 
the solution of © can be calculated by mean of 



P^^{V,t) 



Fa{z)p^{V,t"z)dz. 



(7) 



The formula is especially useful for some particular cases 
whose the exact solutions of the ordinary FPE have a 
closed form (for example, the harmonic potential leading 
to a linear force field in FPE). It is interesting also to 
observe the integral representation of the fractional FPE 
solution in [2^ (see expression (2.32)). Now clearly, that 
formula is nothing else as a consequence of the subordi- 
nation relation ([5]) (or ([7])). In this connection it should 
be noted that the relation ([5]) is not of convolution type, 
so the derivation of Eq. ^ , having the fractional integral 
of time, from the ordinary FPE is not entirely trivial. 

Many papers (23 - [28| have focused on the derivation of 
the fractional FPE with different potentials and its so- 
lution. Starting with (29j . the continuous time random 
walk approach is very popular for that goal. However, 
only recently it has been shown that the solutions are 
density functions of a stochastic process [s^l- We sup- 
port the latter point of view: the problem of anomalous 
diffusion should be analyzed with the exact definition of 
the corresponding random process. The density function 
and the master equation are derived from this process. 



IV. FLUCTUATION-DISSIPATION RELATION 
AND H-THEOREM 

According to Eq. ([l} , the variance of the random vari- 
able V is 

{V^{r))^vloe-'^^ + ^{l-e-'^^), (8) 

where Vi^o is the initial condition. Since the random pro- 
cesses V and S{t) are independent, we average the ex- 
pression ([U on the internal variable t so that 



(V^it))^ / FM){V,'rz))dz, (9) 
Jo 



4 



where the line over a variable denotes the average over 
the internal variable r. Calculating the following integral 



i^„(z)e-2^*°^dz = S„(-27t"), (10) 



the exponential functions in ^ are replaced with the 
Mittag-Leffler function for ([9|). The stationary state of 
(j4|) is finite so that 



lim {V^{t))=D/j. 



(11) 



The constants D and 7 are interpreted as generalized dif- 
fusion and damping coefficients respectively. The mean 
{Vi{t)) is zero as well as (^^(t)) = 0. The boundary case 
a = 1 may be also included in the study, as S(t) = t. 
The probability density p^{t,T) reduces to the Dirac 6- 
function, and Eq. ([T]) becomes the ordinary Langevin 
equation in the true time t. 

In the other hand, the energy of a classical system is 
distributed equally among all degrees of freedom. We 
get the fluctuation-dissipation relation D/j = kBT/m 
for the given temperature of a bath T and the mass of a 
particle m, and ks is the Boltzmann constant. The ex- 
pression is very similar to the Einstein relation, but the 
constants D and 7 are generalized. It should be pointed 
out that the stochastic time arrow does not break this 
equal distribution law and influences only on the char- 
acter of relaxation (slow power decay). Therefore, in 
this case the concept of temperature is valid, i. e. the 
stationary state of the fractional FPE is defined by the 
temperature T. The difference of entropies of equilib- 
rium and arbitrary states gives a Lyapunov functional 
A{t) > 0. No wonder that its temporal evolution con- 
firms the H-theorem. Although the law of relaxation 
toward thermal equilibrium changes, it remains mono- 
tonic, and the equilibrium state has the most entropy 
due to the Gibbs-Boltzmann distribution. The fact, that 
no modifications of the Boltzmann thermodynamics for 
anomalous diffusion described by the equation of type (|6]) 
are required, was already noted in [26, 28, 31]. However, 
the true cause of the result was not established. Now it 
is clear that both processes ([IJ and (|1]) are closely con- 
nected and have a common ground. 



V. GENERAL KINETIC EQUATION WITH 
THE STOCHASTIC TIME CLOCK 



For a general type of a Markovian process the general 
kinetic equation is 



dt 



k=0 



{WnkPk{t)-WknPn{t)} , (12) 



where W^n are the transition probability rates from state 
n to state k. This equation defines the probability p„ for 
the system to be in state n. The term WnkPk describes 



transitions into the state n from states fc, and WknPn 
corresponds to transition out of n into other states k. 
The continuous version of Eq. takes the form 



dP{y,t) 
dt 



{W{y I y')P{y\ t) - W{y' \ y)P{y, t)] dy' 



with the initial condition P{y^O). Let us represent both 
these equations as 



dpjt) 
dt 



Wpit), 



(13) 



where W denotes the transition rate operator. It is im- 
portant to emphasize here that this operator is time- 
independent. The equation (|13[) can be written in the 
integral form 

p{t)^p (0) + f drWp (t) . 
Jo 

The Laplace transform p (s) with respect to t is given by 

/•oo 

p{s)= / e-'*p{t)dt 
Jo 

and leads to 

Wp{s) = sp{s) -p{0). 

Now we determine a new process with the probability 
equal to 

Pa{t) - / /(t,T)p(T)dT. 



In Laplace space the probabilities Pa (t) and p (t) are re- 
lated by Pais) — s"~^p(s"), where 



Pais) 



' Pait)dt 



is the Laplace image of Pait). By the simple algebraic 
transformations we find 



Wpais) = S^-^Wpis") 

= s"-i{s"p(s")-p(0)} 
= s"p,(s)-p(0)s"-i. 

Thus, the fractional extension of Eq. ([T51) reads 



(14) 



Pait) =p(0) 



r(«) Jo 



dTit~T)"-'WpaiT). (15) 



For a = 1 we recover Eq. ([T^ . For a system with discrete 
states the generating function is of the form 



Gic,t) = j2c'pkit), 



k=0 



5 



p^(0)-p,(0) - 
Consequently, we obtain 



1/2 




FIG. 2: The relaxation of the probabilities pf [Ea. (|17p ] and 
Pi [EQ. p8[) ] plotted as a function of time. The dotted line 
corresponds to the equilibrium state. 



(17) 
(18) 



The steady state of the system corresponds to equilib- 
rium, p^(oo) — p^{oo) — 1/2 (Fig. [5]). The transition 
rate w is defined by microscopic properties of the system 
(for instance, from the given Hamiltonian of interaction 
and the Fermi's golden rule). The value (2w)~^/" may 
be interpreted as a generalized relaxation time. The ran- 
domization of time clock essentially changes the character 
of relaxation in such a two-state system. If only a 7^ 1, 
the relaxation has an algebraic decay. In this connec- 
tion it should be mentioned here that the experimental 
relaxation curves of glasses show just the algebraic decay 



where the restriction | C 1^ 1 is imposed to ensure con- 
vergence of the series. With the help of the generating 
function, one can find the moments by taking the deriva- 
tive with respect to C and then setting C = 1- The gener- 
ating function of the process governed by the stochastic 
time clock is given by the relation 



/>oo 

G^iCt) = / F^{z)G{C,ez)dz. 
Jo 



(16) 



Thus, the generating function for a discrete Markov pro- 
cess directed by the process S{t) can be obtained from 
the appropriate generating function of the parent process 
by immediate integration. 

As an example, we consider the relaxation in a two- 
state system. Let TV be the common number of objects 
in this system. If N-^ is the number of objects in the 
state tj is the number of objects in the state I so 
that N = + N^. Assume that for t = the states t 
dominate, i. e. 

N^{t = 0) 
N 



Pt(0) = 1, 



m = 0) 

TV 



where p-^ and p^ are the probabilities to find the system 
in the states t and l respectively. Denote the transition 
rates by w. In the case the general kinetic equation with 
the stochastic time clock (ITSt is written as 



pt(0 = Pt(o) 



piit) = pm 



w 



r(a) Jo 



w 



{t~Tr~'{p^iT)-p^iT)}dr, 



(t-r)"-i{pt(r)-p^(r)}dT 



r(a) Jo 

From the linearity of these equations it follows 



VI. SUMMARY 

We have shown that the fractional FPE can be derived 
by using the ordinary Langevin equation. Although the 
processes described by Eq. ([6]) are non-Markovian at the 
true time, they are Markovian with regard to the internal 
time. So the strange kinetics results in the randomiza- 
tion of time clock of a Markov process. In the proba- 
bility theory the operation is called the subordination of 
one random process by another [l^. As has been stated 
above, the subordination does not break the fiuctuation- 
dissipation relation and the H-theorem. The stochastic 
time clock has a clear physical sense - a particle interacts 
with a bath in random points of time so that there are 
memory effects. It should be noted that the memory is 
a direct consequence of the random time steps belonging 
to the strict domain of attraction of a a-stable distribu- 
tion. One and only one index a characterizes both the 
corresponding a-stable process and its hitting time pro- 
cess. The stochastic differential equation (|4]) describes a 
random velocity field directed by a random Markov pro- 
cess. In this case the dynamical foundation of statistical 
physics is valid. The procedure of the randomization of 
time clock extends the domain of applicability for the 
general kinetic equation. 
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